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ABSTRACT: 

In  this  paper  the  standard  wing  equation,  as  normally  derived  from 
lifting  line  theory,  is  further  refined  and  a  solution  procedure  more 
basic  than  the  usual  collocation  technique  is  developed.  The  calculation 
method  adopted  avoids  the  necessity  of  performing  an  explicit  matrix 
inversion;  all  equations  can  be  solved  sequentially,  one  at  a  time.  On 
the  other  hand  this  technique  involves  the  evaluation  of  numerous  integrals 
over  the  span.  The  calculations  are  cumulative,  and  can  be  carried  as  far 
as  necessary  to  achieve  any  required  degree  of  accuracy.  The  analysis  is 
interesting  not  only  for  purposes  of  practical  calculation  but  also  for 
the  light  it  sheds  on  the  essential  mathematical  structure  of  the  basic 
aerodynamic  phenomena  involved.  This  same  general  method  of  calculation 
can  also  be  readily  adapted  to  the  solution  of  other  common  types  of 
engineering  problems. 
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1.      INTRODUCTION 

The  standard  wing  equation,  as  usually  derived  on  the  basis  of  ele- 
mentary lifting  line  theory,  is  normally  solved  by  the  collocation  technique. 
While  this  approach  gives  acceptable  results  under  normal  conditions, 
circumstances  can  arise  which  make  a  somewhat  deeper  analysis  desirable. 
The  first  purpose  of  this  paper  is  to  show  that  the  standard  wing  equation 
itself  can  with  advantage  be  further  developed  and  refined.    In  particular, 
it  is  profitable  to  split  this  equation  into  two  parts  which  respectively 
govern  the  basic  and  additional  lift  distributions,  and  to  refine  the  result- 
ing pair  of  equations  by  an  integration  process.    The  second  purpose  of 
the  paper  is  to  develop  a  calculation  technique  which  is  somewhat  more 
fundamental  and  powerful  than  the  usual  collocation  procedure.    The  method 
selected  is  based  on  the  theory  of  orthogonal  functions.    It  provides  an 
explicit  solution  procedure  which  can  be  carried  out  as  far  as  necessary 
to  satisfy  the  initial  equations  to  any  required  level  of  accuracy.    The 
solution  obtained  is  equivalent  to  that  obtainable  by  matrix  inversion,  but 
the  present  technique  achieves  this  result  by  an  indirect  process  in  which 
an  explicit  matrix  inversion  is  not  required  at  any  point.    The  development 
of  this  method  is  of  interest  not  only  for  purposes  of  practical  calculation 
but  also  for  the  theoretical  light  it  sheds  on  the  mathematical  structure  of 
wing  aerodynamics.    The  same  general  calculation  technique  can  also  be 
readily  adapted  to  the  calculation  of  other  types  of  engineering  problems 


as  well;  the  wing  problem  represents  only  one  particular  application  of 
the  method. 

All  symbols  are  defined  in  Section  2.     For  readers  interested  primarily 
in  the  resulting  method  of  calculation,  rather  than  in  the  detailed  analysis 
itself,  the  calculation  procedure  is  summarized  in  Section  3.    The  essen- 
tials of  the  standard  wing  equation  are  reviewed  briefly  in  Section  4  and 
the  suggested  refinements  are  presented  in  Section  5.     Solution  by  the 
usual  collocation  method  is  briefly  discussed  in  Section  6.     Orthogonal 
functions,  which  provide  the  basis  for  the  proposed  alternative  method  of 
solution,  are  analyzed  in  Section  7,  and  the  resulting  solution  is  developed 
in  Section  8. 

2.      DEFINITIONS  OF  SYMBOLS 

dCT 
a  wing  lift  curve  slope, L 

da 

a   ,        coefficients  which  occur  in  Eq.   (7.8) 
nk 

a  section  lift  curve  slope 

GO 

a  mean  section  lift  curve  slope  as  defined  in  Eq.   (3.4) 

00 

A  initially  unknown  constants  which  define  the  additional  lift 


n 


distribution,  Eq.   (3.34) 


A*  constants  defined  by  Eq.   (8.11) 
n 

AR  aspect  ratio  2  s/c 

b  elements  in  orthogonality  matrix,  Eq.   (3.12) 

B  initially  unknown  constants  which  define  the  basic  lift  distribution, 

n  Eq.   (3.33) 


c  local  wing  chord  at  spanwise  station  7],  ft. 

c  mean  chord,  S/2s,  ft. 

c  section  lift  coefficient  at  location  71 
I 

C  initially  unknown  constants  which  define  the  circulation  function, 


n 


Eq.    (4.3) 


CT  wing  lift  coefficient 

C  induced  drag  coefficient  of  wing 
i 

f  forcing  function  defined  in  Eqs.   (3.23)  or  (3.24) 

F  forcing  constants  defined  in  Eq.    (3.25) 

I  planform  integrals  defined  in  Eq.   (3.9) 

J  .  wing  integral  defined  in  Eq  (3.11) 
nk 

K  ,  integrals  defined  in  Eq.   (7.6) 

s  wing  semi-span,  ft. 

r  2 

S  wing  area  ,  ft. 


V  remote  velocity ,  ft/sec 

00 

w         downwash  velocity  at  wing  station  71,  ft/sec 

a  angle  of  attack  of  wing,  that  is,  angle  between  zero  lift  direction 

of  mid-span  section  and  V   ,  rad 

<y  angle  of  attack  of  wing  at  which  resultant  wing  lift  vanishes,  rad 

-Li 

O 

a  angle  of  twist,  that  is,  angle  between  zero  lift  direction  at  station 

7]  and  zero  lift  direction  at  mid-span,  a  prescribed  function  of  71 
or  9  ,  rad . 

a  twist  angle  at  wing  tip,  rad 

1 

2  / 
T  circulation  about  wing  at  spanwise  station  71,  ft  /sec 


B  I 

6  )    Induced  drag  factors  defined  in  Eqs.   (3.30),    (3.31)  and  (3.32) 

6 

A 

6  ,  Kronecker  delta ,  defined  in  Eq.    (7.2) 
nk 

T]  spanwise  location  expressed  as  fraction  of  semi-span 

9  angular  coordinate  defined  in  Eq.   (3.1) 

S.     basic  and  additional  lift  functions  defined  in  Eqs.   (5.25)  and  (5.26) 

A  orthogonal  lift  function  defined  in  Eq.    (8.23) 
n 

v  planform  function  defined  in  Eq.   (3.6) 

a  auxiliary  parameter  defined  in  Eq.   (3.5) 

t  twist  function  defined  in  Eq.   (3.7) 

t  mean  twist  defined  in  Eq.   (3.8) 

t  relative  angle  of  zero  lift  as  defined  in  Eq.   (3.27) 


o 


n 


orthogonal  functions  defined  in  Eqs.   (7.2)  and  (7.4). 


Y  wing  functions  defined  in  Eq.   (3.10) 

n 


Y  '         wing  functions  defined  in  Eq.   (4.14) 


[      basic  and  additional  downwash  functions  defined  in  Eqs.   (5.27) 
and  (5.28) 

"B  J 

Subscripts 

i,  k,  n      indices  which  may  take  on  the  values  3,5,7 N 

A  pertains  to  additional  lift 

B  pertains  to  basic  lift 


Matrix  Symbols 

f    ) 

V         y  column  vector  containing  elements  3,5,7 N 

I    J 

*\         V  column  vector  containing  elements  3,5,7 n  <N 

v        /  n 

j  r 

>         row  vector  containing  elements  3,  5,  7,  --  N 


a  square  matrix  with  rows  and  columns  corresponding  to 
indices  i  =  3,  5,   7, N. 


A  square  matrix  with  rows  and  columns  corresponding  to 

indicesi  =  3,5,7, n<N. 

n 

the  identity  matrix 


3.      SUMMARY  OF  THE  CALCULATION  METHOD 

The  proposed  method  of  calculation  is  summarized  in  full  in  this 
section.    Only  such  explanations  are  included  here  as  are  essential  for 
making  the  method  of  calculation  clear.    All  questions  of  derivation  and 
interpretation  are  deferred  to  other  sections  of  this  paper.    The  formulas 
are  listed  in  the  approximate  order  in  which  they  would  be  used  in  the 
calculations. 
Preliminaries 
Tj  =  cos  9         Change  of  variable  (3.1) 

—         g 

c  =    —  Average  chord  (3 . 2) 


2s 

AR  =    — 


Aspect  ratio 


(3.3) 


(I 


a     =  i  a 

CO  00 


If) 


dT|       Mean  slope  of  section  lift  curve 


(3.4) 


CT   =    a°° 


7TAR 


An  auxiliary  parameter 


(3.5) 


v(9)  = 


T(9) 


/  a  °°\   /  c 


ale 


a 


.1 


°t 
ll' 


t     =    1    TudTl 
o 

-7T/2 


Planform  function 


Twist  function 


Mean  twist 


(3.6) 


(3.7) 


(3.8) 


I      = 
n 


u  sinnBd  9  n  =  3,  5,   7,  —  N       Planform  integrals     (3.9) 


4      .        -  ,  /  sin  n  9 

Y     =   -  smn9+nau      — ; — -  -   I 
n        7T  1     sin  6         n 


n  =  3,  5,  7,  — -  N 


Wing  functions 


(3.10) 


ff/2 


'nk 


Y    Y    d< 
n    k 


n,  k  =  3,  5,  7  —  N      Wing  integrals      (3.11) 


J 


Orthogonality  Matrix 

The  matrix  [b]  is  calculated  as  summarized  below.    It  has  the  structure 


[b] 


33 


53 


73 


B 


55 


N3 


75 


N5 


77 


N7 


bNN; 


(3.12) 


The  first  three  rows  are  calculated  as  follows 


Row  n  =  3 


-2 


33 


33 


(3.13) 


Row  n    =  5 


a53   "   "b33  J53 


(3.14) 


-2  2 

b55      "    J55  "  353 


b53   "  b55  353   b33 


(3.15) 
(3.16) 


Row  n   =  7 


73 


75 


h   "2     =    T 
77  J77 


33 


53  55 


73 
f75 


73 


=   b 


Fb33 
^53 


55 


(3.17) 


(3.18) 


(3.19) 


Any  desired  additional  number  of  rows  can  be  calculated  by  repeated 
application  of  the  following  sequence,  namely 
Row  n 


n3      1 


n5 


< 


V 


n,n-2i 


b33         ° 


b53         b55 


'  b      _    ^  b      _    ^ 
.     n-2,3     n-2,5 


J 


n3 


j    Jn5 


bn-2,n-2        Jn,  n-2y 

j  ^ 


(3.20) 


n3 


n3 


b"2  =  j 

nn  nn 


;    a  a   _ 

;      n/  ^     n7 


V    n,n-2J    ^  n,n-2 


(3.21) 


fb 


t      r 


n3 


n3 


n5     I     ,  -     n5 

L=  b       s 
nn  ^ 


Vbn/n-2J 


n,n-2^ 


33 


53 


j     n-2,  3 


55 


n~2,5 


0 
0 

b 


-\ 


(3.22) 


n-2,n-2| 


Basic  and  Additional  Lift 


To  calculate  the  additional  lift  distribution,  set 


f  (6)    =   v  (9)  -  -  sin  6 


Additional  forcing  function  (3.23) 


To  calculate  the  basic  lift  distribution,  set 

f  (9)  =  u  (6)  [t  (6)  -  t  ]  Basic  forcing  function  (3.24) 

and  change  the  constants  A    below  to  B    . 

n  n 


f/2 


F     =    !  f(9)  t     (9)  d9 
n         I  n 


J 


o 


n  =  3,  5,   7,  —  N 


Basic  or  additional  forcing  (3.25) 

constants 


Solution  Constants 


fA    1  f  F     • 

^  r3 


A. 


K 


T  5 

►  =  [b]     [b]     <    F7    ► 


AW 


!    — 

I 

*   F 


J 


(3.26) 


Wing  Characteristics 

N 

°Lo 


nB  I         -    t 
n  n 


Zero  lift  angle 


(3.27) 


n   =  3 ,  5 ,  7 , 


N 


a    / 


1=1  + 


•  (  1  +I  "Vn  ) 


Wing  lift  curve  slope 


n  =  3,  5,  7 


(3.28) 


D. 


1    \  r 

7TAR  L      °°  t] 


)2    6B  +  (a.VL)6BA+C'    <1  +  8a) 


(3.29) 


Induced  Drag 


where 


N 


6     = 

B 


nB 


n 


n=3,  5,  7,  — 


k  = 


N 

Induced 

V      nB  A 
;           n  n 
L 

r 

drag 

n=3  ,  5  ,   7 ,   - 

i 

factors . 

N 

2 

i 

i 

6     = 

A 


nA 


n 


n=3,  5,   7,— 


Lift  Distribution 


N 


IB    j_c_\ 

a"  a        — 
"tjl  c  j 


77 


B      sin  n 
n 


n=  3,   5,   7,~ 


Basic 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


JLIc 


CLlc 


77 


sin  6  + 


N 


A      sin  Tl  9 
n 

n=3 ,5,7,-- 


Additional 


(3.34) 


4.      REVIEW  OF  THE  STANDARD  WING  EQUATION 

By  utilizing  the  standard  lifting  line  idealization  of  a  wing,  and  em- 
ploying the  Biot  Savart  law,  we  can  deduce  [1 ,  2]  *  that  the  downwash  w 
at  any  point  T]  along  the  span  is  related  to  the  distribution  of  circulation 
r along  the  span  as  follows: 


*  Numbers  indicate  references  listed  at  end  of  paper. 
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+1 

d  /    r 


dTl' 


w         1  diyuv  s, 

v"  -  ;  l    "  ' —  (4-1} 


J       en'  -  n) 


-i 

It  is  customary  to  introduce  the  change  of  variable 

T|  =  cos  e 

(4.2) 
T,'  =    cos  91 

and  to  describe  the  distribution  of  the  circulation  Tby  a  series  of  the  form 


/    r 

4V   s 


CO 


=    Y  cn  sin  n6'  (4-3) 


n  =  1,  3,  5,  — 

The  present  discussion  is  restricted  to  lift  distributions  which  are 
symmetrical  about  the  mid- span  and  hence  only  terms  of  odd  index  are 
retained  in  Eq.    (4,3). 

Upon  introducing  Eq.   (4.3)  into  Eq.   (4.1)  and  making  use  of  the 
integral  formula 


1.  cos  n  9'   d  8'  sin  n  9 

77  J  .,  sin  9 

cos  9'    -  cos  9 


the  following  solution  is  obtained  for  the  downwash 


(4.4) 


w 
V    , 

00 

V      r              n  sin  n  9 
L         n              sin  9 

n=l,3/5/ 

(4.5) 


11 


Once  the  above  integration  has  been  performed,  there  is  no  further 
need  to  distinguish  between  the  variables  T]  and  T]1  or  9  and  9'  and  the 
prime  marks  may  be  deleted  from  Eqs.    (4.2)  and  (4.3). 

The  lift  exerted  by  an  element  of  the  wing  lying  between  T]and  (T|  -*-  dT|) 
may  be  expressed  firstly  in  terms  of  the  section  lift  coefficient  c,  and 
secondly  in  terms  of  the  circulation  F,  in  accordance  with  the  Kutta 
Joukowski  law.     Equating  these  two  alternative  expressions  for  the  local 
lift  gives 


c 


/£     MS(^I  (4.6) 

I  C;  <       co      ' 


where 


—     =    average  chord  (4 .  7) 

<is 


2s 
£R    =   —   =    aspect  ratio  (4.8) 


The  local  lift  coefficient  c«   can  be  expressed  in  terms  of  the  local 
section  lift  curve  slope  a     and  the  local  net  angle  of  attack,  taking  into 
account  the  effect  of  the  local  downwash  angle    —  \  in  reducing  the  effec- 

00  ' 

tive  angle  of  attack.      Hence  we  obtain 


W       X 

2.      =     a       {a    +    a     -  — -)  (4.9) 


'I  co        V  t  V 


00 


where  a  =  wing  angle  of  attack  =  angle  of  attack  of  the  zero  lift 
direction  of  the  airfoil  section  at  midspan  with  respect 
to  the  remote  velocity  V    . 

CO 

12 


and  at  (9)   =  angle  of  twist  =  angle  of  zero  lift  direction  at  T]with  respect 

to  zero  lift  direction  at  midspan,a  prescribed  function  of  0. 
Upon  combining  Eqs.   (4.3),    (4.5),    (4.6)  and  (4.9)  and  rearranging 
we  can  finally  obtain  the  fundamental  wing  equation  in  the  following  form, 
namely, 

)         C     V   (9)    =   av  {ot  +   a)  (4.10) 

L,  n   n  t 

n  =  1,3,5, 

where 


c 
a  a  t'ldT]    =     mean  section  lift  curve  slope  (4.11) 


fi 


'C 

! 
u  0 


a    =       °°       =     an  auxiliary  parameter  (4.12) 


7TAR 


v  (9)      =[  z~H"3"\    =     a  prescribed  planform  function  (4.13) 


a     l  \   c 

00 


.    /« \  4  „  n  sin  n  9  , .    ,  „v 

Y'    (9)    =  —     sin  n  P    +    au    — ~. — : —  4.14 

n  7T  sin  fi 


Apart  from  minor  changes  in  notation,  Eq.   (4.10)  is  the  fundamental 
wing  equation  as  found  in  most  elementary  texts  on  aerodynamics  [1,2]. 
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5.      ALTERNATIVE  EQUATIONS  FOR  BASIC  AND  ADDITIONAL  LIFT 

By  substituting  Eq.   (4.3)  into  (4.6)  multiplying  through  by  dT]    = 
-  sin  6  d  0,  integrating  over  the  semi-span  and  taking  advantage  of  the 
orthogonality  relation 

7T/2 

j 

sin  n  9  sin  k  9  d  9  =  6  ,       =     +   1  if  n  =  k 

nk 

J  0 

=    0  if  n^  k  (5.1) 

we  obtain  a  solution  for  the  overall  lift  of  the  wing  in  the  form 
1 

CL    =     1        CJ[-)dTl  =    7rmCi  (5'2) 

Jo 

This  verifies  the  well  known  fact  that  the  overall  lift  of  the  wing  is 
fixed  solely  by  the  first  term  of  the  series  in  Eq.   (4.3) .     The  other  terms 
affect  the  form  of  the  lift  distribution  but  not  the  resultant  net  lift. 

It  is  useful  to  express  the  twist  distribution  of  the  wing  in  dimension- 
less  form  as  follows 


a 
_t 

<y. 


t  (9)  =  —      =    a  prescribed  twist  function  (5.3) 

V 


where    a      =    angle  of  twist  at  wing  tip  (5.4) 

1 

The  basic  wing  relation,  Eq.   (4.14)  ,  can  be  further  developed  by  in- 
troducing the  substitution 


14 


c 

n 


a   a      B      +  CT    A 
00  t       n  L      n 


(5.5) 


7T  ,AR 

n   =  1,  3,  5,  ■ 


The  constants  B     now  define  the  basic  lift  distribution  of  the  wing, 
n 

This  is  the  distribution  associated  with  the  wing  twist.     Specifically  it 

is  the  distribution  which  occurs  when  the  wing  is  at  the  angle  of  attack 

<y  =   <y        which  gives  C     =  0.    The  constants    An    define  the  additional 

o 
lift  distribution  that  occurs  as  a  result  of  changes  in  angle  of  attack. 

Hence  it  is  proportional  to  the  angle  {at  -   o/T    )  or  to  CT  itself. 

Lo  L 

Consider  Eq.    (5.5)  for  the  specific  case  n  =  1.     Substitute  for  C. 
from  Eq.   (5.2).    The  resulting  expression  must  be  satisfied  for  all  values 
of  CT  .    In  particular  by  setting  C,  =  0  we  find  that 

Bl  =  ° 

A    =  1  (5.6) 


Substituting  the  expressions  (5.5)  into  the  basic  wing  relation  Eq. 
(4.10),  and  making  use  of  the  special  result  (5.6),  we  obtain 


a<y       Y      B     r    (6)  +-^|   sine+YA     r  (P) 
t     /        n    n  ttARL  /      n    n 

n  =  3,5,7, n=3,5,7, 


(Jv  (ot+  a)  (5.7) 


This  must  be  satisfied  for  all  operating  conditions.     In  particular  for 

a  =  a         and  C     =0,  it  reduces  to 
L0  L 
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a    )      B    Y'     (0)      =    av  (a       +  a)  (5.8) 

t,  L        n    n  L  t 

n  =  3 ,5 , 7, 


a 
t 
1 


Then  subtracting  Eq.    (5.8)  from  (5.7)  we  obtain 


7TAR1 


I    sinG+y   A  i'   (9)    I    =    au(ff-ct    )  (5.9) 

n  n         _l  LQ 

n  =  3  ,  5  ,  7, 


Eq.    (5.8)  can  be  further  reduced  by  dividing  through  by  o<y      and 


1 


setting 


^f)    =    TL0  (5-10> 

'  1 

Eq.   (5.9)  can  be  further  reduced  by  substituting 

0L-a(.-«L)  (5.11) 


/dC 
where     a    =f L\      =     wing  lift  curve  slope  (5.12) 

da 


In  this  way  we  obtain  for  the  basic  and  additional  lift  the  two  funda- 
mental equations. 


V      B     Y     (9)     =     p(tt       +    t)  (5.13) 

2,       n    n  LQ 

n  =  3,5,  7, _ 

sin  8+y  A     y'   (9)    =      -^    u  (5.14) 

L      n     n  \    a    / 

n  =  3 ,5 ,7 , 
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The  complete  analysis  of  a  twisted  wing  can  be  obtained  either  by 
solving  the  standard  wing  relation,  Eq.    (4.10)  for  two  different  angles  of 
attack,  or  alternatively,  by  solving  the  above  pair  of  relations,  Eqs.   (5.13) 
and  (5.14)  ,  for  basic  and  additional  lift.     The  total  number  of  unknowns  to 
be  found  is  exactly  the  same  by  both  methods;  only  the  form  is  different. 
It  can  be  said,  however,  that  Eqs. (5. 13)  and  (5.14)  are  somewhat  more 
fundamental  in  form  in  the  sense  that,  for  a  wing  of  prescribed  geometry, 
the  unknowns  B    ,   r n.  A     and!     °°\  are  true  invariants,  whereas  each  one 

n    °   n     \r\ 

of  the  "constants"  C     in  Eq.    (4.10)  is  actually  a  function  of  angle  of 

n 

attack. 

Moreover,  Eqs.   (5.13)  and  (5.14)  lend  themselves  to  further  develop- 
ment in  a  useful  and  interesting  way.    In  this  connection,  the  first  step 
is  to  multiply  these  equations  through  by  dT|  =    -sin  9d9  and  integrate 
over  the  semispan. 

This  operation  involves  the  following  integrals,  Firstly, 

1  1 

( 
'  a 

00  \  I  c 

udT]    = 


■J 


a    /  \  c 

00 


dTl  =    1  (5.15) 


0  0 

a  result  which  follows  directly  from  the  definition  of    a     as  given  in 

GO 

Eq.    (4.11). 

ly, 

T7/2 

sin  n9  sin9d9  -  1  for  n  =  1 

=  0  for  n  /=  1       (5.i6) 
0 
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Secondly, 

1 

r 

4 

sin  nQ  d7?    = 

4 

77 

Thirdly,  we  define  the  average  value    t    of  the  twist  function  by  the 

integral 

1  tt/2 

( 


T    =      \    TudTj  =       I        to  gin  6d9  (5.17) 

J  J 

0  0 

Finally,  we  define  the  series  of  integrals 
1  tt/2 

I      =  u    Sm  "  6     dT\  =  (5.18 

n  i         sin  6 


n  =  1,3,5, 


We  also  note  from  this  last  definition  that  specifically  for  n  =  1 
1 


l1    =  \     udT]  =    1  (5.19) 


0 

Upon  integrating  Eqs.    (5.13)  and  (5.14)  over  the  semi-span,  utilizing 
the  foregoing  integrals,  then  simplifying,  we  obtain  the  two  significant 
results  below,  namely, 


t         =     >     nB  I      -    t  (5.20) 

LQ 

n  =  3,5,7,  
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a 

00 

=     1    +    Q       1 

a 


V  nA  I 
/.        n  n 


(5.21) 
n  =  3,5,7, 

These  results  are  useful  for  two  reasons.     First  of  all,  they  can  be 

used  to  eliminate  the  unknowns    t        andf_°f\    from  Eqs.   (5.13)  and  (5.14) 

0  la    ) 

thus  reducing  the  number  of  simultaneous  unknowns  which  must  be  found 

from  these  equations.     Secondly,  they  provide  a  subsequent  solution  for 

these  two  unknowns,  and  one  that  is  especially  accurate  because  it  is 

based  on  an  integration  over  the  entire  wing.  _ 

/3oo     \ 

After  Eqs.    (5.20)  and  (5.21)  have  been  used  to  eliminate  tt    and; 1 

LQ         \    a    / 

from  Eqs.    (5.13)  and  (5.14),  and  after  the  results  have  been  simplified, 
they  become  finally 

/  Bf    (8)    =    d(t-  t)  (5.22) 

n  n 

n  =  3,5,7, 

V        A   i|r      (9)     =    (v  -  £    sin  o)  (5.23) 

)  n   n  tt 

n  =  3,5,7, 

where 

t    (9)  =   ;  sin  n9+  naul  **in  "9    -    I    )  (5.24) 

n  Tt  x  sin  9  n/ 

n  =  3,5,7, 

Eqs.   (5.22)  and  (5.23)  represent  an  alternative  form  of  the  funda- 
mental wing  equations.     They  can  advantageously  replace  the  standard 


19 


wing  equation  (4.10)  as  usually  employed  in  the  lifting  line  analysis. 
The  series  on  the  left  side  of  Eqs.    (4.10)  and  (5.22)  and  (5.23)  are,  of 
course,  truncated  at  some  value  n   =  N,  depending  on  the  accuracy 
required.     Often  only  the  four  terms  in  C.  ,  C    ,  C     and  C     are  retained  in 

-LOO  / 

Eq.   (4.10).       Note  that  in  this  case,  Eqs.   (5.22)  and  (5.23)  will  contain 
only  three  unknowns  apiece,  namely,  B    ,  B    ,  B     and  A    ,  A    .  A    .     More- 

0  0/  00/ 

over  the  latter  solutions  can  be  expected  to  be  somewhat  more  accurate 

since  they  include  the  averaging  effects  of  integration  over  the  semi-span, 

Note  that  Eqs.    (5.22)  and    (5.23)  are  of  identical  format  on  the  left. 

Only  the  forcing  functions  on  the  right  are  different.     Eq.    (5.22)  shows 

clearly  that  the  constants  B    depend  only  on  the  dimensionless  planform 

and  twist  functions  and,  in  particular,  that  an  untwisted  wing  gives 

B     =0.     Likewise,  Eq.   (5.23)  shows  that  the  constants    A      depend  only 

on  the  deviation  of  the  planform  function    v    from  an  ellipse  of  equal  area 

and,  in  particular,  that  a  wing  having  an  elliptical  planform  function 

gives  A    =0. 
n 

It  should  also  be  noted  that  while  Eqs.   (5.22)  and  (5.23)  are  more 
effective  than  Eq.   (4.10)  in  the  sense  noted  above,  they  entail  the  addi- 
tional labor  of  evaluating  the  integrals  Tand  I 

n. 

After  Eqs.   (5.22)  and  (5.23)  have  been  solved  for  the  constants  B 

*L  \ 
and  A    ,  the  angle  of  zero  lift  /       0  \  =    t        and  the  slope  of  the  wing 

1 


lift  curve     —      can  then  be  found  accurately  from  Eqs.   (5.20)  and  (5.21), 
respectively. 
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The  resulting  basic  and  additional  lift  and  downwash  distributions  can 
be  found  by  combining  Eqs.   (4.3),   (4.5),   (4.6)  and  (5.5).    We  thereby 
obtain  the  following  solutions. 


/c 


IB     \ 


a   ot 

00  t. 


S.\ 


i  ■  c 


-  y      B     sin  n  fi   =    \     (p) 
ff  _j        n  B 

n  =  3 ,5  ,7, 


[5.25) 


XA 


c  i 


7TAR 


a    a 

°°  t. 


B       _     V 


w    \ 


7TAR 


77 


sin 


+  /    A     sin  n  6       |  = 
Li     n 

n  =  3,5,7,-- 


V9) 


n 


n  sm  n 
sin  9 


°fe 


(6) 


n  =  3,5,7, 


(w    \ 
A  i  V      a  n  sin  n  6 

7^—  1  =     1+)     A        — : — - —      =    uja  a 
V    /  Li        n  sin  6  A 

CO  ' 


n  =  3,5,7  , 


(5.26) 


(5.27) 


(5.28) 


Once  the  constants  A    and  B    are  known,  they  fix  the  four  functions 

n  n 

above  which  define  the  fundamental  aerodynamics  of  the  particular  wing 

under  consideration. 

The  one  remaining  question  of  basic  significance  which  falls  within 

the  purview  of  lifting  line  theory  concerns  the  distribution  and  overall 

value  of  the  induced  drag.    The  overall  coefficient  of  induced  drag  is 

defined  by  the  integral 

1  )■ 


Di 


«m^  = 


0 


c  c 

c  A     ~  +   c  *       ~ 

i_  -  I.    ~ 

c  Ac 


w           w. 
B             A 

V       +    V 

00                              CO 

dH 

(5.29) 
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Upon  substituting  for  the  quantities  in  Eq.    (5.29)  in  terms  of  the 

four  functions  defined  in  Eqs.    (5.25)  through  (5.28),  then  rearranging,  we 

! 

i 

obtain 

7T/2 

2 


a   at     \  \ 

CT  \   °°     V  A_  (G)  m    (p)     sin  0d9 

Di  7TAR 


0 
77/2 


a   a     CT 

00  t,      L 


+ 


k—      I  F*«  <e>  »,  <•>  +  X.  («  %  (fl 


sin  9d6 


77  AR  Mb'    a  w       av/b 

1 

J 

0 

77/2 


c2 

+    ItAR         I      XA(9)  V^         Sined"  (5'30) 

0 


The  integrals  on  the  right  side  of  Eq.   (5.30)  show  that  there  are  three 
distinct  components  of  induced  drag  and  the  respective  integrands  define 
the  corresponding  distributions.    However,  instead  of  a  direct  evaluation 
of  these  three  separate  integrals,  it  is  simpler  to  proceed  as  follows. 
First  we  evaluate  the  overall  integral  in  Eq.   (5.29)  using  Eqs.   (4.3),    (4.5) 
and  (4.6).    The  result  is  simply 

v      2 

n        =    77AR   )      nC  (5.31) 

D.  '  n 

i 

n  =  1,3,5  , 
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Next  we  substitute  Eq.   (5.5)  into  (5.31),  expand  and  regroup  terms.     In 
this  way  we  obtain  finally 


c  l 


d.  it  ar 

i 


/-  2 

la.*t   CT.^na    +    C 


6b  +  \^\clYba  +  cl  (i  +  v_;(5-32) 


2 

where        I        =  Y      nB,_  (5.33) 

B  [_  u 


6 
BA 


n  =  3,5,7, 

=     Y    nB  A  (5.34) 

/  n  n 


6 
A 


n  =  3,5,7, 

2 


=    V     nA  (5.35) 

/  n 

n  =  3,5,7, 


We  also  note  in  passing  that  for  the  special  case  of  an  aerodynamically 
untwisted  wing  of  elliptical  planform,  Eq.   (5.26),   (5.28),   (5.21),  and 
(5.32)  reduce,  respectively,  to  the  following  familiar  and  elementary 
forms,  namely, 

c    l-\  4 

I  \—     =    CT    —  sin  6  =    elliptical  lift  distribution  (5.36) 

I  C   /  L     77 

CL 
—      =     — —     =     uniform  downwash  along  span  (5.37) 


a 


a      = 


i  +  a< 


slope  of  wing  lift  curve  (a  maximum)  (5.38) 


77  AR 
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C  =      — rr-      =     induced  drag  (a  minimum)  (5.39) 

i 

These  last  four  results  illustrate  the  fact  that  the  shift  from  the 

standard  form  (4.10)  of  the  wing  equation  to  the  revised  forms  (5.22)  and 

(5.23)  does  not  affect  the  essential  character  of  the  solution  obtained. 

6.      THE  COLLOCATION  METHOD 

The  usual  method  of  obtaining  an  approximate  solution  of  the  standard 
wing  relation,  Eq.    (4.10)  is  by  collocation,  that  is,  by  satisfying  this 
equation  exactly  only  at  a  limited  number  of  discrete  points  along  the 
semi-span.     The  number  of  points  must  of  course  equal  the  number  of 
unknown  constants  C,  ,  C„,  C    , C       to  be  retained  in  the  approximate 

J.  o  o  IN 

solution.    The  location  of  the  collocation  points  is  somewhat  arbitrary, 
but  it  is  customary  to  space  them  at  equal  intervals  of  9  along  the  semi- 
span,  to  include  the  midspan  point  at  9=— ,  and  to  omit  the  wingtip  point 
9=0.     The  solution  of  the  resulting  set  of  simultaneous  equations  entails 
a  matrix  inversion,  but  is  otherwise  quite  straightforward.    If  a  computer 
be  utilized,  of  course,  even  the  matrix  inversion  poses  no  appreciable 
computational  burden. 

In  practice  the  above  calculations  even  if  performed  only  on  a  desk 
calculator,  usually  work  out  in  a  reasonably  satisfactory  manner.    This 
is  partly  owing  to  the  fact  that  any  well  designed  wing  of  moderate  to  high 
aspect  ratio  has  a  lift  distribution  that  does  not  deviate  strongly  from  the 
optimum  elliptical  distribution.     Consequently  the  successive  constants 
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C    ,  C    ,  C    , typically  diminish  rapidly  in  magnitude,  and  only  a  few 

terms,  often  just  four,  suffice  to  produce  results  of  an  accuracy  adequate 
for  many  engineering  purposes.     Moreover,  a  quest  for  extreme  numerical 
accuracy  would  not  be  realistic  because  the  basic  lifting  line  idealization 
on  which  the  wing  equation  is  based  is  itself  an  approximation. 

Nevertheless,  there  are  situations  in  which  the  mathematical  circum- 
stances are  not  so  favorable,  for  example,  in  cases  of  strong  discontinuities 
in  the  lift  distribution  owing,  say,  to  the  action  of  wing  flaps  over  part  of 
the  span.    Moreover,  apart  from  practical  questions  of  computation,  col- 
location has  a  number  of  features  which  are  not  optimal  from  the  theoretical 
viewpoint.    For  example,  the  choice  of  the  actual  collocation  points  in- 
troduces an  element  of  arbitrariness  into  the  solution.    It  would  be  in 
some  respects  preferable  to  have  a  method  whose  accuracy  depends  in  a 
known  way  only  on  the  number  of  terms  retained  in  the  solution,  and  not 
on  such  contingencies  as  the  location  of  collocation  points.    Furthermore, 
the  collocation  method  can  under  certain  adverse  circumstances  lead  to  a 
set  of  equations  which  are  numerically  ill  conditioned.    Moreover,  the 
purely  theoretical  conception  of  how  the  collocation  solution  approaches 
its  limit  as  the  number  of  terms  retained  in  the  solution  hypothetically 
increases  without  bound  tends  to  be  somewhat  obscure;  it  entails  the 
conceptual  difficulty  of  visualizing  the  inversion  of  a  matrix  in  the  limit 
as  the  number  of  rows  and  columns  increases  without  bound. 
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7.      ORTHOGONAL  FUNCTIONS 

There  is  an  alternative  to  the  collocation  method  that  deals  with  the 

above  difficulties,  both  conceptual  and  computational,  in  a  very  clear  and 

decisive  way.    It  can  be  developed  through  the  introduction  of  a  certain 

family  of  orthogonal  functions  $     (9),  n  =  3,5,7, that  are  related  in 

n 

a  simple  way  to  the  basic  functions  if    (9)  that  occur  in  the  wing  equation. 

The  mathematical  structure  of  these  auxiliary  functions  is  developed  in 

this  section.     It  turns  out  that  the  functions  $    ,  while  essential  to  the 

n 

theoretical  development,  can  be  eliminated  from  the  final  calculations; 

only  certain  constants  that  occur  in  the  definition  of  these  functions  carry 

over  into  the  resulting  calculation  procedure. 

The  auxiliary  functions  $     (9)  are  related  to  the  original  functions 

n 

if    (9)  by  equations  of  the  following  form: 
n 

V9'    "    b33  *3 


'5(6)    =    b53*3    +    b55+5  (7.1a) 


(9)    =    b       ilr      +    b       i+b       i 
7  K  '  73  ¥3  75  T5  77  v7 


>    (9)    =    b   _  *     +  b   _  if.    + +b       if  (7.1b) 

n  n3     3         n5    5  nn    n 


The  constant  coefficients  b  ,    in  these  equations  are  chosen  so  that 

nk 

the  functions  $    (9)  form  an  ortho-normal  set,  that  is,  so  that 
n 
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7T/2 


i     $.    d6  =    t>  .     =    0       for  k    =    3,5,7,  — -  (n-2) 
n    k  nk 


=    +1     for  k  =  n 


(7.2a) 
(7.2b) 


It  is  convenient  to  develop  much  of  the  theory  in  matrix  format.    Ac- 
cordingly, we  introduce  the  following  abbreviated  matrix  notation. 
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0 

[b]n  = 

b73 

b75 

11 
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0 
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n7 

-- 

b 

nn 

_ 

Consequently,  Eqs.    (7.1)  may  be  written  simply  as 
$ 


=  tb]n  ♦ 

•  n  I  Jn 


(7.3) 


(7.4) 


where  we  call  [b]  the  orthogonality  matrix. 

In  connection  with  the  ensuing  analysis,  it  will  also  be  convenient 
to  define  the  following  auxiliary  integrals,  namely, 
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and 


tt/2 
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K 


nk 


fn  (6)  $k(9)  dP  k  =  3,5,7,  —  (n-2) 


(7.6) 


The  integrals  K       will  ultimately  be  eliminated  from  the  analysis. 
The  integrals  J  ,,  however,  can  be  assembled  to  form  the  following  matrix, 
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(7.7) 


It  can  be  seen  that  the  matrix  [b]      „  fixes  the  functions  <L  ,  $_  ,  $,,, 

n-z  3       5       7 

$      0.     Now  suppose  that  this  matrix  is  known  and  that  we  wish  to  add 

n-z 

to  it  the  initially  unknown  n      row  of  coefficients  which  will  define  the 

next  function  $   .     For  this  purpose,  it  is  useful  to  write  $    temporarily  in 
n  .  n 

the  alternative  form 


n-2 


>       =    b 
n  nn 


a         $       +  t 
nk      k  n 


k  =  3,5,7, 
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(7.8) 


In  matrix  notation,  this  may  be  written  in  either  of  the  two  equivalent 
forms  below,  namely, 
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Also,  from  (7.4) 
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Now  premultiplying  Eq.   (7.10)  through 


by  ^  if         gives 
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(7.13) 


We  now  integrate  this  over  the  semi-span,  that  is,  from  9  =  0  to 


TT 


2,  then  apply  the  ortho-normal  condition  (7.2)  and  the  definitiions  (7„5) 
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and  (7.6).     The  following  simplifications  occur,  namely, 
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(7.15) 


where  [I]    _9  is  the  identity  matrix,  with  rows  and  columns  corresponding 
to  indices  3,  5,   7, (n-2). 
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With  the  aid  of  relations  (7.14)   (7.15)  and  (7.16),  the  integrated 
form  of  Eq.   (7.13)  can  be  rearranged  to  yield  the  following  significant  result 
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This  defines  the  solution  for  the  constants  a  .    such  as  to  satisfy  the 

nk 

orthogonality  conditions  (7.2a) .     However,  a  further  reduction  of  the  right 
side  of  Eq.   (7.17)  is  possible  and  useful.     Note  that,  from  the  definition 
(7.6)  and  from  (7.4), 
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Combining  (7.17)  and  (7.18)  gives  the  orthogonality  relation  in  the 
desired  final  form,  namely, 

\   }  (1 

<%>        =-[bln-2]'J  (7-19) 


v      'n-2  l  n-2 


To  apply  the  normalizing  conditions  (7.2b),  we  first  post-multiply 

-2 
(7.9)  by  (7.10),  then  multiply  through  by  b    ■  .      This  gives 
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n  n  I  i 
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We  next  integrate  this  over  the  semi-span  noting  that,  from  (7.2b), 

7T/2 

e2  d9   =    1  (7.21) 
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and  from  (7  .  5) 
tt/2 

r 
i 

i  2 

i       tlr    de 

J    n 


=  J. 


nn 


(7.22) 


Also,  the  first  term  on  the  right  side  of  (7.20)  may  be  simplified  after 
integration  in  the  manner  of  (7.15).     The  next  two  terms,  which  incidentally 
are  of  equal  magnitude,  can  be  simplified  according  to  (7.16)  .    In  this 
way  we  ultimately  obtain  from  the  integrated  form  of  (7.20)  the  result 


I    1 


b  =    \a 
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This  reduces  readily  to 
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(7.24) 


We  can  eliminate  <  K    r  from  this  by  using  the  orthogonality  rela- 

•      n-  2 
tion  in  the  form  given  in  (7.17).     The  final  result  therefore  becomes  simply 
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=     J 
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1     n 


(7.25) 


n-2 


At  this  point  the  three  relations  (7.10),   (7.19)  and  (7.25)  together 
define  the  function  $    in  such  a  way  that  it  satisfies  the  normalizing  con- 
dition (7.2b)  and  is  orthogonal  to  all  the  previous  functions  in  the  set, 

namely,  to  f_,   $_  ,  $  _ , $     n.    However,  it  still  remains  to  determine 

3      5      7  n-z 

the  coefficients  b   _,  b   r,  b   _ b  „  which  are  needed  to  describe 

n3       n5       n7  n,  n-2 

$    in  the  standard  format  of  Eq.    (7.4) .    In  this  connection  we  note  first 
n 

that 
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n5  f 


n7 


V    nn, 
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nn         n 
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Now  a  comparison  of  (7.26)  with  (7.9)  shows  that 


(7.26) 


-T      f 
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=     b  a 
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(7.27) 


Moreover,  by  analogy  with  Eq.   (7.4)  we  may  write 

f 


^     n-2 


=  [b]n     9    I     tin     9  (7'28) 

n-2   V.     >  n-2 


Substitution  of  (7.28)  into  (7.27)  discloses  finally  that 
T  T 


J>  b     f      =     b       \  a     V  W      0  (7.29) 

n-2 


n-2  v         n-2 


This  is  the  required  result  which,  along  with  (7.19)  and  (7.25), 

J.L. 

enables  the  n       row  of  matrix  [b]     to  be  calculated  when  all  the  previous 

n 

rows  are  known. 

The  principal  results  of  this  section  can  now  be  summarized.    It  is 

supposed  that  all  quantities  over  the  range  of  the  index  i  =  3,5,7, 

(n-2)  are  known  from  previous  calculation.    The  following  calculations 

then  extend  these  results  to  the  additional  index    i  =  n. 

Firstly,  we  must  evaluate  the  additional  integrals 
7T/2 


Jnk 


(6)     t    (6)dP  k  =  3,5,7,— n  (7.29) 

n  k 


J 

0 


For  any  given  wing  planform,  the  functions  ^  as  defined  by  Eqs. 

(5.24)  are  known.    Hence  the  constants  J  ,    are  computable  by  direct 

nk 


numerical  integration. 
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Next  we  calculate,  in  the  order  given,  the  quantities 


a  =-  [h]  ;j     I  (7.30) 

n  n-2  >   n   r 

\ 

/ 
n-2  n-2 

,  T        , 

-2 
b        =    J         -      a  a  (7.31) 

nn  nn  n  n 

n-2  n-2 


,T  /        ]T 


'  b  =    b      J  a  [b]  (7.32) 

n  nn    ■    n  t  n-2 

n-2  n-2 

Eqs.   (7.31)  and  (7.32)  fix  row  n  of  the  orthogonality  matrix  [b]    . 

n 

Consequently,  the  function  <$    is  now  defined  in  the  standard  format,  namely, 

n 


$        ='    b  S  Y  (7.33) 

n  n 

n  n 

Moreover,  the  function  $     so  defined  satisfies  the  ortho-normal  rela- 

n 

tion  (7.2)  with  respect  to  all  previous  functions  in  the  set. 

The  above  calculation  sequence,  Eq.   (7.29)  through  (7.32),  may  then 
be  repeated  for  the  next  successive  row  of  the  matrix  [b] .     In  this  way, 
the  fundamental  matrix  [b]  can  be  expanded  to  any  desired  number  of  rows 
and  columns . 

We  shall  next  show  how,  once  the  orthogonality  matrix  has  been  es- 
tablished to  any  desired  size,  all  remaining  unknowns  of  the  problem  can 
be  calculated  to  the  corresponding  degree  of  accuracy  in  a  routine  way. 
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8 .      SOLUTION  BY  MEANS  OF  THE  ORTHOGONALITY  MATRIX 


Both  of  the  revised  wing  relations,  Eqs.   (5.22)  and  (5.23),  are  solved 
by  the  same  general  method  and  hence  for  the  purposes  of  this  section,  it 
is  preferable  to  adopt  a  generalized  notation  which  embraces  both  equations 
Also,  it  is  helpful  to  retain  the  matrix  notation.     Hence  we  rewrite  Eqs. 
(5  ..22)  and/or  (5.23)  simply  as 
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Ar 
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(8.1) 


8.2) 
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For  calculating  the  additional  lift  distribution  we  set 


f(0)    =    fA(9)      =     u    -    ~  sin 


(8.3) 


For  calculating  the  basic  lift  distribution  we  set 


f(e)    =   fB(e)     =     v    (t-  t) 


(8.4) 


In  this  latter  case  we  also  change  ^  A^  to  j  Bpor  consistency  with  the 
notation  established  earlier. 


It  is  also  convenient  to  define  the  set  of  integral  constants 
7T/2 


F      = 
n 


f(9)  f     (P)  dP 
n 


n  =  3,5,7,  ---  N 


(8.5) 
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and,  of  course,  the  corresponding  vector 


3 

! 


Also  since 

f     I 


V 


it  follows  that 


4=  [b]"M  $f 


and  that 


T  T 
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-T 


and  that 


T  T 

1     \  _T    \ 

=  i  $  V      [b]    '  <    A 


Mr 


t 


(8.6) 


(8.7) 


(8.8) 


(8.9) 


(8.10) 


This  last  result  suggests  the  usefulness  of  defining  the  auxiliary 


vector 

<   A*  /•      =         [b] 


!aI 


(8.11) 


j     -  n 

We  now  find  that  the  basic  relation  to  be  solved,  Eq.   (8.1)  ,  can  be 


rewritten  in  the  alternative  form 


J 


(8.12) 


•(  AH      =        f(9) 

I  J 

The  great  advantage  of  (8.12)  over  the  original  form  (8.1)  lies  in  the 
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orthogonality  properties  of  the  functions  $   .    For  if  we  now  multiply  (8.12) 

through  by  any  one  of  the  functions  $  (6) ,  then  integrate  over  the  semi- 

n 

span,  we  find  that,  because  of  the  orthogonality  properties,  all  integrals 

but  the  one  in  n  vanish  on  the  left  with  the  simple  result  that 

tt/2 
v 

i 

A*       =  f(9)  $     (P)d9  (8.13) 

n  n 

0 

n    =    3,5,7, (n-2) 


Upon  substituting  into  (8.13)  the  relation 
,      .T  -      , 


-        i 
*n    =       MR 

L     jnt 


(8.14) 


and  utilizing  the  definitions  (8.5) ,  this  result  can  be  rewritten  in  a  form 
more  convenient  for  calculation,  namely, 


■ 


A*     =    <  b      f     \    F  f  (8.15) 

n 


n 


This  can  readily  be  generalized  over  the  entire  vector  <  A*>  .       The 
process  amounts  to  pre-multiplying  Eg.   (8.12)  through  by  (8.7),  then  in- 
tegrating over  the  semi-span  and  applying  the  ortho-normal  conditions  in 

the  usual  way.    We  thereby  obtain 

■•  (      > 

i  i       I 

\  A*  ■       =       [b]\  F   ^  (8.16) 
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A  significant  property  of  the  constants  A*    as  revealed  by  Eqs.    (8.15) 

n 

or  (8.16)  is  that  they  are  true  invariants  of  the  solution  in  the  following 

sense.     The  value  of  A*     depends  on  F.  ,  F_,F_, F    .     But  owing  to  the 

n  6       5      /  n 

zeros  in  [b] ,  it  is  independent  of  all  the  higher  order  terms  F         n  ,     F 

n  +  2         n  +  4  , 

and  so  on.     Consequently,  adding  on  additional  terms  to  the  solution  in  no 

wise  affects  the  values  of  the  constants  A*     previously  calculated. 

n 
(     * 
Nevertheless,  while  -j  A    "  contains  the  true  invariants  of  the  solution, 

j 

it  is  more  convenient  in  some  ways  to  continue  to  express  the  solution  in 

terms  of  the  original  unknown  constants,  as  contained  in  the  vector   '  A  >. 

( 

To  solve  for    {  A  £     we  need  merely  to  premultiply  Eq.   (8.11)  through  by 

T 

[b]     .     The  result  is  simply 

j  A  ■   -    [b]T  ]  A*  I  (8.17) 

In  many  cases  it  may  be  preferred  to  bypass  the  vector  i  h* t  altogether 
in  the  solution.    To  show  how  this  can  be  done,  we   simply  substitute  for 
from  (8.16)  into  (8.17).    The  result  becomes,  finally, 


{a*} 


I-) 


ik}=    [b]T[b]{] 

Thus,  once  the  orthogonality  matrix  [b]  is  known,  Eqs.   (8.5)  and  (8.18) 
summarize  the  solution  for  the  initially  unknown  constants  A   .    When  these 
are  fixed,  all  remaining  features  of  the  solution  can  be  established  from 
the  equations  developed  in  the  earlier  sections. 

The  solution  represented  by  Eq.   (8.18),  although  it  entails  no  explicit 
matrix  inversion  at  any  point,  is  in  fact  equivalent  to  a  matrix  inversion. 
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To  show  this,  we  premultiply  the  basic  relation,   Eq.   (8.1)  ,  through  by 
<^£  then  integrate  over  the  semi-span  and  apply  the  definitions  of 
J       and  F     as  given  in  Eqs.   (7.5)  and  (8.5)  respectively.     In  this  way 


we  obtain 


[J]  4A>     =   JF  >  (8.19) 


The  solution  of  this  equation  is 

i      1  -iJ       1 

Af=      [J]  F   f 

Comparison  of  Eq.   (8.18)  with  (8.20)  discloses  that 


(8.20) 


[b]T[b]      =      [J]"1  (8.21) 

which  demonstrates  the  sense  in  which  calculation  of  matrix  [b]  amounts 

to  the  inversion  of  matrix  [J]  . 

Bearing  in  mind  that  the  elements  of  [b]  are  calculated  in  the  first  place 

solely  from  the  elements  of  [J] ,   it  is  clear  that  what  we  have  here  in  effect 

is  a  particular  method  of  inverting  matrix  [J].     But  since  [J]  could  be  any 

symmetrical  matrix,  the  method  amounts  to  an  inversion  procedure  applicable 

to  any  symmetrical  matrix.    While  the  method  happened  to  be  developed 

in  this  instance  in  terms  of  the  orthogonal  functions    $    ,  these  functions 

n 

do  not  actually  occur  in  the  resulting  calculation  formulas.     Hence  the 
final  algorithm  may  be  applied  to  invert  any  symmetrical  matrix  quite  ir- 
respective of  its  origins.    In  this  particular  instance,  the  elements  of  [J] 
happen  to  be  certain  specific  integrals;    the  same  inversion  procedure 


40 


should  also  apply  to  symmetrical  matrices  whose  elements  are  defined  in 
some  other  way.    Hence  the  technique  is  clearly  applicable  not  just  to  the 
wing  problem,  but  to  a  wide  range  of  other  applications  as  well. 

One  advantage  of  solving  by  means  of  the  matrix  [b]  is  that  the  cal- 
culations are  strictly  progressive  and  cumulative.     The  process  of  adding 
on  each  additional  row  to  [b]  provides  additional  finer  detail  about  the 
solution  without  disrupting  the  details  calculated  from  earlier  rows.     This 
becomes  somewhat  clearer  if  the  solution  be  expressed  in  terms  of  the 

constants  A*  rather  than  in  terms  of  the  original  constants  A    .     The  in- 
n  n 

variance  properties  of  the  A*  have  been  noted  earlier.    In  contrast  with 

n 

this,  the  constants  A     fluctuate  each  time  a  higher  order  term  is  added 

n 

to  the  solution,  and  approach  their  true  ultimate  values  only  in  the  limit 
as  the  number  of  terms  increases  without  bound.    This  analysis  therefore 
sheds  much  light  on  the  nature  of  the  convergence  process  that,  in  prin- 
ciple, defines  the  exact  mathematical  solution  of  the  original  basic 
equation. 

The  solution  may  therefore  be  viewed  as  the  result  of  superimposing 
an  indefinite  number  of  modes  each  of  which  is  orthogonal  to  all  of  the 
preceding  modes.     To  illustrate  this  concept,  consider  the  solution  for 
the  additional  lift  distribution,  Eq.    (5.26),  rewritten  in  matrix  notation, 
that  is , 


T 
A    /  -c  \         4  .         4 


— -    l-\=     -    sin    9   +   -\  A  >    \    sin  n9  > 
°L     lc]         *  f  j(  J 


I 

(8.22a) 
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rr 
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(8.22b) 
(8.22c) 
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5        i 
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1   Xn(9)>    H  V  (9)V=    ff     M    N     Sln  ? 
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(6)  j 

N      y 
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(8.23) 


sin  N  e 


J 


Thus  the  functions  ^    (9)  ,   X    (9)  ,  X  (9) , X    (9)    as  defined  in  Eq. 

(8.23)  represent  the  successive  invariant  orthogonal  modes,  a  linear 
combination  of  which  fixes  the  additional  lift  distribution  as  shown  in 
Eq.   (8.22c) .     The  addition  of  any  mode  n  to  the  solution  does  not  affect 
the  amplitudes  of  the  modes  3,  5,  7,  --  (n-2)  previously  obtained.     This 
example  illustrates  the  theoretical  concept,  at  any  rate.    Of  course,  for 
practical  calculation,  the  ordinary  form  shown  in  Eq.   (8.22a)  would 
probably  be  preferred  as  involving  somewhat  less  overall  computation. 
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